Draft version February 1, 2008 

Preprint typeset using I£T^X style cmulatcapj v. 7/15/03 



o 
o 

(N 

Or 
<D 



6 

CO 



> 

in 
o 

ON 

o 
l> 
o 



THE SURVIVAL RATE OF EJECTED TERRESTRIAL PLANETS WITH MOONS 
John H. Debes 1 ,Steinn Sigurdsson 2 

Draft version February 1, 2008 

ABSTRACT 

During planet formation, a gas giant will interact with smaller protoplanets that stray within its 
sphere of gravitational influence. We investigate the outcome of interactions between gas giants and 
terrestrial-sized protoplanets with lunar-sized companions. An interaction between a giant planet and 
a protoplanet binary may have one of several consequences, including the delivery of volatiles to the 
inner system, the capture of retrograde moons by the giant planet, and the ejection of one or both 
of the protoplanets. We show that an interesting fraction of terrestrial-sized planets with lunar sized 
companions will likely be ejected into interstellar space with the companion bound to the planet. The 
companion provides an additional source of heating for the planet from tidal dissipation of orbital and 
spin angular momentum. This heat flux typically is larger than the current radiogenic heating of the 
Earth for up to the first few hundred million years of evolution. In combination with an atmosphere 
of sufficient thickness and composition, the heating can provide the conditions necesary for liquid 
water to persist on the surface of the terrestrial mass planet, making it a potential site for life. We 
also determine the possibility for directly detecting such systems through all-sky infrared surveys or 
microlensing surveys. Microlensing surveys in particular will directly measure the frequency of this 
phenomenon. 

Subject headings: planets and satellites:formation-astrobiology-methods:N-body simulations 



1. INTRODUCTION 

Binary protoplanets of terrestrial mass and lower may 
be common-large Kuiper belt objects in particular seem 
to have a high fraction of satellites that formed through 
giant im pacts much in the s ame way as the Earth/Moon 
system (|Canupl l200ll2005| ). Three of the four known 
large KBOs (Pluto, 2003UB313, and 2003EL61) contain 
satel lites, with Pluto/C haron having a mass ratio closest 
to 1 (|Brown et alJ l2~606). Assuming similar densities and 
albedos between parent and body and based on the rela- 
tive brightnesses, the other two binary KB Os have mass 
ratio s similar to the Earth/Moon system ([Brown et al.l 
2006). Since these objects and the Earth/Moon system 
exist, three body interactions need to be studied in the 
context of planet formation. 

Giant impacts, the main process for forming lunar 
mass moons to terrestrial planets, occur 1-2 times over 
the course of terrestrial planet formation simulations 
with enou gh angular momentum to form th e Earth-Moon 
system ([Agnor et al.l 119991 : IChamberdl200l . The value 
of this angular momentum is dependent on the angle of 
impact, the mass of the impactor, and ratio of impact 
speed to escape speed. A range of parameters can cre- 
ate a situation where Earth-Moon analogues are formed. 
This range is limited for the specific case of the Moon 
due to oth er constraints such as a lack of iron relative to 
the Earth ([Jones fc Delanolll989l ). but many constraints 
are relaxed for an impact that creates a generic moon of 
lunar mass, allowing for a reasonable probability that the 
Earth/Moon system is mediocre and analogues should be 
present in a non-negligible fraction of planetary systems. 

Even though there arc examples of binarity in the ter- 
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restrial planet region and in the Kuiper belt, it is not 
clear whether conditions exist to form Earth/Moon sized 
systems where giant planets arise. Two scenarios could 
create a situation where a binary protoplanet may come 
in close contact with a Jovian planet. 

The first is that the binary interacts with a Jo- 
vian in situ during giant planet formation. The bi- 
nary would be a leftover from oligarghic growth, simi- 
lar to the scenario proposed by the formation of Uranus 
and Neptune through scat tering with Jupiter and Sat- 
urn (Th ommcs et~aLl[2002l ). In some cases, oligarchic 
losers will beejected to escape speeds. The size distri- 
bution of protoplanets in this region is poorly known, 
but oligarchic growth starts at ~ 10~ 5 M e and pro- 
ceeds to_j L _gas_acn^ion_mnaw^v_gro^ mass of ~10 
M ffi (|Thommes et al.l 120031: iPollack et all 1 19961) . Oli- 
garchic growth in the terrestrial planet region appears 
to spread from highe r density to lower density regions 
(|Kokubo fe Idal I2002T > . which suggests that a size range 
of oligarchs would be present in the region of 5-10 AU. 
Leftovers that range from the smallest oligar ch mass to 
ice giant mass are possible ([Chiang et al.ll2007h and could 
explain Uranus's peculiar axis of rotation as a result of 
a gia nt impact between it and a large late stag e left- 
over ((Kor^cansk^eEaD QUO; HTHninLeraD [2002) . The 
constraints on the size of left-overs at late stages of the 
Solar System's formation around the region of Neptune 
places stringent constraints on the maximum size of the 
planetesimal distribution to sizes much smaller than the 
Earth, though one time e ncounters with Mars- like proto - 
planets are not ruled out (|Murrav-Clav fe ChianeJl2006l) . 
However, the typical architecture of a planetary system is 
not well known and the stochastic behavior of oligarchic 
growth would argue for a wide range of possible outcomes 
for planetary systems between 5-10 AU, including the 
prese nce of terrestrial mass protoplanets ([Levison et al.l 
1998). 
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Secondly, forming planetary systems where one or 
more giant planets migrate to small orbital radii may 
lead to the ejection of sever al smaller planets in the pro- 
cess ([Raymond et al.ll2006[ ). Some ejected planets may 
have bound companions that are a significant fraction of 
their host planet's mass. The frequency of such systems 
and the survival rate of these binary protoplanets during 
ejection is an interesting open question. 

If terrestrial planets are ejected with lunar-sized com- 
panions, they may be habitable, despite being deprived 
of an insolating central star. If an isolated terrestrial 
planet is ejected with significant atmosphere captured 
from the parent protoplanetary disk, the planet will have 
a su rface temperatu re high enough to sustain liquid wa- 
ter (|Stevensonlli999ft . 

The presence of the companion will add heat through 
tidal dissipation of spin and angular momentum, which 
will augment the heat flux due to radioactive materials 
in the interior of the planet. Tidal heating has been 
investigated in the melting of the moon, the vulcanism 
on Io, and the heating of Europa dPeale fc Cassenlll978t 
iPeale et all [19791: iCassen et all fl979F iCarr et all 11998ft . 
The heating will decay over time, but if life evolves, it 
could p ersist for long times below the surface of the larger 
planet (|Laughlin fc Adams! 12000ft . 

2. NUMERICAL SIMULATIONS 

To test the frequency of binary protoplanet ejections 
we have run 2700 num erical simulatio ns using the Mer- 
cury integration code (| Chambers! fl 999ft . We placed a bi- 
nary protoplanetary system with an Earth mass primary 
and a lunar mass secondary in Hill unstable orbits rel- 
ative to a giant planet on a concentric orbit. Hill sta- 
bility is defined as a system that experiences no close 
appr oaches between the consit uent planets in the sys- 
tem (jHilll 118861 : lGladmanlll993ft . We chose initial posi- 
tions for the terrestrial planet to be between 0.5-1 times 
the critical Hill radius interior and exterior to the giant 
planet. Initial eccentricities of the binary protoplanet 
system were set to 0.05 and inclinations were < 5°. Ini- 
tial semi-major axes of the companion's orbit were set 
to 30 R®, half the current Moon-Earth separation. Ini- 
tial Moon-Earth separations would be smaller, but we 
chose this separation as a compromise between simula- 
tion speed and physical reality, as this separation dic- 
tates the time step required to accurately integrate the 
equations of motion. Timesteps were set to ~8 days and 
were adaptable. A system was run for a maximum of 
10 6 yr but was checked at 10 4 yr and 10 v r. In this way, 
highly interacting systems that disrupted the binary pro- 
toplanets were removed early in the integration. Our cri- 
teria for ejection was that the binary protoplanets were 
bound to each other, but unbound from the central star 
with a radial distance > 15 AU from the central star. 
In the simulations, no system was considered ejected un- 
less it passed beyond a radial distance of 5000 AU, so 
in most cases the ejected binary protoplanets were al- 
lowed to evolve for several thousand years beyond the 
perturbation that ejected the system. Fractional errors 
in energy and angular momentum were ~10 -11 but were 
larger if systems suffered a collision. The highest errors 
were ~10 -8 . Figure Q] shows a typical ejection of a pro- 
toplanet with a surviving companion. 

As a further test of the stability of the simulations we 



ran ten long term simulations of an isolated binary pro- 
toplanet system to determine if the initial orbit could be 
stable for 10 6 yr, the maximum length of the simulation. 
We found that in all ten simulations the orbits remained 
stable for initial eccentricities of ~0.15. 

3. RESULTS 

Of the 2700 simulations, 90, or 3.3% of the interactions 
resulted in only the Earth being ejected, while 123 or 
4.6% of the simulations ended in an ejection of a bound 
Earth-Moon system. In most cases the eccentricity of 
the bound Earth/Moon systems increased after ejection. 
Bound ejections occur through multiple grazing encoun- 
ters with the Jovian, which is consistent with the pertur- 
bations on the Earth/Moon system being less than their 
binding energy. 

An estimation of the tidal dissipation in these ejected 
systems shows that tidal heating can be as important as 
radioactive heating. Heating can come from circulariza- 
tion of a satellite's eccentric orbit as well as synchroniza- 
tion between the planet's spin and the satellite's mean 
motion. Heating rates for circularization are dependent 
on the eccentricity e, semi-major axis a, and mean mo- 
tion n: 
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The variable fi is the ratio of elastic to gravitational 
forces, or the effective rigidity of the planet. The vari- 
able Q p is the planet's specific dissipation function, and 
R p is planetary radius. Heating for synchonization can 
be found by the following equation: 
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where we assume that the planet's rotational frequency 
is greater than its mean motion and t he rotational period 
for t he terrestrial planet is 10 hr (|Murrav fc Dermottl 
2000) . Median values of eccentricity and semi-major axis 
of the ejected pairs are 0.13 and 30 R®, respectively. 

We assume two cases, a rocky planet and an icy planet. 
A protoplanet will fall in between these two cases, de- 
pending on where it formed. A rocky planet will be iden- 
tical to the Earth in mass, radius, density, and rigidity. 
We t ake fc 2 =0.299, /i=4 and Q p =12 (jMurrav fc Dermottl 
2000). For an icy body we assume a rigidity of ice of 
4x 12 dyne cm~ 2 , and assuming a density of 1 g cm~ 3 
for pure water ice, a /c2p=0.7 and /2=1. We assume a 
Q P =100. 

From among the end states of the ejected planets, the 
maximum heating came from a system with a semi-major 
axis of 21 R ffi and an eccentricity of 0.21, corresponding 
to a total heating of 4.2xl0 22 ergs/s assuming a rocky 
planet. This is 10 times large r than current terrestrial 
radiogenic heating (IStacevlll992ft . The maximum heating 
for an icy planet is a factor of five larger. The distribu- 
tions of heating rates for rocky and icy planets from our 
simulations relative to the terrestrial radiogenic rate are 
presented in Figures [H 

We compare these rates to th e tidal heating on Io, 
whose luminosity is ~10 21 ergs/s (|Veeder et al.lll994ft . or 
about three times larger than the terrestrial radiogenic 
heating rate. The highest heating of a rocky planet is 
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roughly 4% of the total insolation the young Earth ex- 
perienced from the Sun, assuming a current rate of 1.5 
xlO 24 ergs s^ 1 . While this heating is small in absolute 
terms compared to the conditions present on Earth for 
the formation of life, it does not create an insurmount- 
able barrier for life to form and evolve. The heating 
would provide energy for maintaining liquid water and 
forming life on the planet, potentially in localized hot 
spots where vulcanism or other geothermal processes oc- 
cur. 

Heating decay over time for the terrestrial planet would 
depend on the change in the planet's rate of rotation 
(to) and the semi-major axis (d) due to tidal cvolu- 
tion of the pair. Assum ing d oc a~ 5 ' 5 and ui oc a~ 6 
(|Murrav fe Dermottl [2000h . we calculated a and w as a 
function of time and determined when E in Equation [2] 
meets current radiogenic levels. For our best case, we 
found that the heating from Equation [2] reaches radio- 
genic heating in 144 Myr for a rocky planet, and 246 
Myr for an icy planet. This may be enough time for 
life to arise and adapt to decreasing temperatures on the 
planet. 

4. DETECTABILITY OF FREE FLOATING EARTHS 

Given the uncertainties in the formation rate of binary 
protoplanets, the number of such objects per star, and 
the number of planetary systems, we rely on an analogue 
of the Drake Equation to estimate the space density of 
ejected terrestrial planets and ejected binary protoplan- 
ets. In reality there will be a size distribution of such 
objects each with different rates, but this will allow us 
to get an idea of detectability through imaging and mi- 
crolensing surveys. 

N A r * j 'planet s^planets j 'binary j 'ejected (3) 

where N* is the number of stars in the Galaxy, f p i an ets is 
the frequency of planetary systems around stars, n p i anets 
is the number of terrestrial sized protoplanets that are 
formed in the system, ibinary is the frequency of pro- 
toplanets that have companions, and i e jected is the fre- 
quency these protoplanets are ejected with their compan- 
ions intact. With the exception of N* and iejected, these 
quantities are not known, but can be estimated within an 
order of magnitude; so a reasonable estimate is possible. 

We assume that the fraction of stars that have planets 
is 50% as a compromise between the minimum of a few 
percent, which would be the fraction of solar type stars 
with known planets, and a maximum of one planetary 
system per star. We assume that the number of terres- 
trial protoplanets available for ejection is on the order of 
10, which is the same order of magnitude of rocky mate- 
rial present in the giant planet region. We assume that 
there are several contenders to become the eventual large 
giant planet, and the "losers" form the primary resevoir 
of unstable protoplanets that will encounter the Jovian. 
Systems with Jovians that migrate will also play a role, 
which may help to increase this resevoir. 

The fraction of protoplanets that have companions is 
also unknown, but can be assumed to be roughly 33% 
based on our own Earth-Moon system and the assump- 
tion that our Solar System is mediocre. Considering the 
frequency of binary asteroids and kuiper belt objects, the 
frequency is probably no less than a few percent. 



With these assumptions, we arrive at a total number 
of free floating binary planets of 7 x10 s , assuming a total 
of 10 10 stars in the galaxy. Assuming that these planets 
are restricted to the galactic thin disk, with a radius of 
25 kpc and a scale height of 0.4 kpc, we estimate a space 
density for these objects of 9xl0~ 4 /pc 3 . 

To estimate the free floating planets' spectrum, we as- 
sume a median heating luminosity of 5.3x 10 21 ergs/s for 
rocky planets and a Planck spectrum. The flux density 
will peak at ~77/im, so one would require an all-sky IR 
survey sensitive to ~0.3 nano-Jy to detect ~100 objects 
out to a distance of 50 pc. The background limited sen- 
sitivity of Hershel's PACS instrument is predicted to de- 
tect 3mJy sources at ~68//m. A similar sensitivity would 
require a > 750m size primary to detect a free floating 
planet in 1 hour of observing. However, the closest of 
these objects will be 5 times closer and thus 25 times 
brighter with a large proper motion. Provided it can be 
resolved from the IR background due to galaxies, one to a 
few detections could be made with less stringent restric- 
tions. Thus, we consider the above to be a conservative 
estimate of what could be detected. 

In the case of microlensing, proposed next generation 
sur veys will be sensitive to free floating terrestrial plan- 
ets (|Bennett fc Rhic 2002|). In most cases, microlensing 
events of planets in orbit around stars have a degener- 
acy between mass, distance, and proper motion. This 
degeneracy is fortunately broken given certain observa- 
tional set ups which allow the o ther parameters to be 
measured (|Han et all 12004. 2005). Binary protoplanets 
with separations similar to those studied in this work 
should be easily detectable given that the orbits reside 
well within the Einste in ring radius of the host planet 
([Bennett fc Rhid 12002). Simulations of detection fre- 
quencies for next generation missions assuming a free 
floating planet frequency of 1-10 free floaters per star 
result in ~ 10-100 detec tions of earth mass fre e floating 
planets for the mission (Bennett & Rhic 2002). Assum- 
ing a third of free floating planets started out with com- 
panions and ~5% of those were ejected with a companion 
intact, there would be a 1-2% chance that any free floater 
would have a lunar mass companion. 

5. OTHER IMPLICATIONS 

There are other implications that deserve further 
study. A significant fraction of stable binary systems at 
the end of our simulations were at large radii (~1%), 
but had not been ejected. They will either continue 
to be on wide orbits, or they will eventually be ejected 
on timescales longer than 1 Myr. If some of these sys- 
tems remain bound to the central star, they may become 
rogue planets that periodically perturb planetesimals at 
larger orbital separations, akin to the "Planet X" hy- 
pothes is presented to explain per iodic extinctions on the 
Earth ( Whitm ire fc Mat ese 1985). For binary protoplan- 
ets, this occurrence appears to be rare. It will be more 
common for singular protoplanets, but we did not follow 
this specific outcome. 

Of the 2700 simulations 12% of the lunar companions 
were captured into resonant retrograde orbits with the 
Jovian planet. Tidal dissipation between the compan- 
ion and the Jovian could be enough to capture the com- 
panion into a stable bound orbit. Three body interac- 
tions have been pointed to as a formation mechanism 



4 



Debes et al. 



for the retrograde Neptunian satellite Triton, and could 
explain how the Saturnian satellite Phoebe could have 
been captured from the outer solar s ystem into its cur- 
rent elliptical and retrograde orbit (|Johnson fc Luninel 
[20051: lAgnor fc Hamilton! l200l . 

Several simulations had remnant protoplanets that 
were kicked into eccentric orbits that strayed into the 
inner system, with periastrons that approached 1.5 AU, 
thus providing a mechanism to bring volatile rich proto- 
planets into the inner system. The delivery of volatiles 
from further out in the proto-planetary disk seems to be 
a natural outcome of plane t formation (jMorbidelli et alJ 



U987tlSchar3l2006f) . These objects will also be detectable 
through microlensing. 

Since all of the ejected binary protoplanets experience 
heating greater than that currently experienced by the 
Earth, they will require t hinne r atmo spheres than thos e 
considered bv lStevenso nl (|1999h . From lStevensonl (|1999f ). 
the surface temperature of a protoplanet with an atmo- 
sphere is 

/ r \ 0.36 

425x A ( ) (4) 



0.001 



planet 
. 2004) 



2000; iRavmond et al 

Tidal heating is greater around planets that are more 
massive, especially terrestrial companions to gas giants 
or brown dwarfs. These objects will either be large oli- 
garchic losers or isolate d objects that forme d with large 
circumstellar disks (i.e. iLuhman et aT1l2005[ ). The heat- 
ing these objects feel will be more substantial than that 
felt for the planets in this study and could be longer 
lasting if multiple companions are formed whose mutual 
gravitational interactions provide eccentricity perturba- 
tions. Each of these larger planets will have respective 
habitable zones where tidal heating repla ces insolation as 
the dominant heat source for liquid water (Revnold s~et al.l 



where x is the ratio of tidal heating to the current ra- 
diogenic heating of the Earth, f a tm is the fraction of the 
planetary mass that is in the form of an atmosphere. To 
maintain a surface temperature of 270 K with radiogenic 
and median tidal heating together requires two times less 
atmosphere by mass than radiogenic heating alone. Since 
tidal heating decays rapidly, radiogenic heating will be 
an important mechanism for sustaining an established 
biosphere for longer time scales. 

The authors thank J. Chambers for helpful discussions 
on oligarchic growth and the referee, D. Stevenson for 
suggestions on improving the paper. 
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Fig. 1. — (Top) An integration of a binary protoplanet in a Hill unstable orbit with a gas giant. The thick blue curve is the gas giant, while 
the dashed red curve represents the binary protoplanets. The thin dashed lines represent the periastron and apastron of the binary pair. 
(Bottom) The same integration but in the rest frame of the terrestrial mass member of the binary protoplanets, showing the semi-major 
axis of its companion as a function of time. The eccentricity of the companion remains roughly constant at ~0.15 
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Fig. 2. — The distribution of tidal heating for the 123 ejected binary protoplanets. The solid line curve corresponds to the assumption 
of rocky planets, while the dashed line corresponds to the assumption of icy planets. The vertical line is the amount of radiogenic heating 
currently on Earth. 



